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ABSTRACT 

We describe narrowband and spectroscopic searches for emission-line star forming galaxies in the 
redshift range 3-6 with the 10 m Keck II Telescope. These searches yield a substantial population of 
objects with only a single strong (equivalent width 3> 100 A) emission line, lying in the 4000 - 8500 A 
range. Spectra of the objects found in narrowband-selected samples at A ~ 5390 A and ~ 6741 A show 
that these very high equivalent width emission lines are generally redshifted Lya A 1216 A at z ~ 3.4 
and 4.5. The density of these emitters above the 5ct detection limit of 1.5 x 10~ 17 ergs cm~ 2 s _1 
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roughly 15,000/D°/unit z at both z ~ 3.4 and 4.5. A complementary deeper (1 a ~ 10 
s _1 ) slit spectroscopic search covering a wide redshift range but a more limited spatial area (200 □") 
shows such objects can be found over the redshift range z = 3 — 6, with the currently highest redshift 
detected being at z = 5.64. The Lya flux distribution can be used to estimate a minimum star formation 
rate in the absence of reddening of roughly 0.01 M Mpc -3 yr -1 (Hq — 65 km s _1 Mpc^ 1 , qo = 0.5). 
Corrections for reddening are likely to be no larger than a factor of two, since observed equivalent widths 
are close to the maximum values obtainable from ionization by a massive star population. Within the 
still significant uncertainties, the star formation rate from the Lya-selected sample is comparable to that 
of the color-break-selected samples at z ~ 3, but may represent an increasing fraction of the total rates 
at higher redshifts. This higher- z population can be readily studied with large ground-based telescopes. 

Subject headings: cosmology: observations — early universe — galaxies: evolution — galaxies: formation 



1. INTRODUCTION 

The search for high-redshift galaxies and the effort to 
map the star formation history of galaxies have progressed 
rapidly in the last several years as magnitude-limited spec- 
troscopic surveys pushed into the z = 1 — 3 range (Cowie 
et al. 1996; Cohen et al. 1996), while color-based selec- 
tion techniques produced many objects in the z — 2.5 — 5 
range (Steidel et al. 1996a, 1996b; Franx et al. 1997; Stei- 
del et al. 1998), particularly in the exquisite Hubble Deep 
Field (HDF) sample (Lowenthal et al. 1997). However, the 
galaxies chosen by these techniques correspond to objects 
with ongoing massive star formation and small amounts 
of extinction, and may represent only part of the popu- 
lations at these early epochs. More evolved objects may 
be heavily dust-reddened and more easily picked out at 
submillimeter wavelengths (Smail, Ivison, & Blaine 1997; 
Hauscr et al. 1998), while earlier stages in evolution may 
have relatively little continuum light and be too faint to 
be seen in the magnitude-limited samples or selected with 
the color-break techniques at current sensitivity limits. 

This latter class of objects may represent the earliest 
stages of the galaxy formation process, in which substan- 
tial amounts of metals have yet to form. These galaxies 
may have much stronger Lya emission relative to the stel- 
lar continuum, since they have massive star formation that 



can excite the Lya emission line, but without so much dust 
that the line is suppressed, and this can result in very high 
observed equivalent widths in the range of 100-200 (1+z) 
A (e.g., Chariot & Fall 1993). Such objects may be hard 
to pick out with color-break techniques but be detectable 
in Lya searches of sufficient depth (Cowie 1988; Thommes 
1996). An increased incidence of strong Lya emission does, 
indeed, appear in the color-break samples at fainter con- 
tinuum magnitudes (Steidel et al. 1998). 

Earlier blank-field Lya surveys (e.g., Thompson, Djor- 
govski, & Trauger 1995; Thompson & Djorgovski 1995) 
failed to find such objects, as Pritchet (1994) has summa- 
rized, but their sensitivity lay at the margin of where such 
objects would be expected in significant numbers (Cowie 
1988; Thommes & Meisenheimer 1995). However, Hu & 
McMahon (1996), using very deep targeted narrowband 
searches (la — 1.5 x 10~ 17 ergs cm~ 2 s _1 ), found z ~ 4.55 
Lya-emitting galaxies with the very strong emission and 
weak or undetected continuua predicted for early star- 
forming objects. The advent of 10 m telescopes, along 
with this successful detection of Lya emitters, prompted 
us to undertake a new survey that has been successful 
in detecting blank-field high-z Lya emitters. The present 
Letter describes the early results of this search, which uses 
extremely deep narrowband filter exposures taken with 
LRIS (Oke et al. 1995) on the Keck II telescope to search 
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for emission-line populations at extremely faint flux levels 
(1 a = 3 x 10~ 18 ergs cm -2 s _1 ). This survey picks out 
sources of extremely high equivalent width (W\ > 100 A) 
emission lines as candidates, and then uses followup LRIS 
spectroscopy (§2) to determine if these are Lya emitters. 

The first results from the Hawaii survey with a 5390/77 
A filter, corresponding to Lya emission at z ~ 3.4, were 
described in Cowie & Hu 1998 (hereafter, Paper I), and 
yielded a number of candidate high-z galaxies similar to 
the redshift z ~ 4.55 Lya-emitting galaxies found by Hu & 
McMahon (1996) and emitters at z ~ 2.4 (Pascarelle et al. 
1996; Francis et al. 1997) in targeted searches. In Paper I, 
we showed that the use of color selection on emission-line 
selected objects of high equivalent width (> 100 A) picked 
out objects with continuum colors similar to those of color- 
selected Lyman break galaxies with measured redshifts of 
z ~ 3.4. They also recovered the one field object whose 
previously measured redshift placed Lya within the filter 
bandpass. However, the emission-line galaxies selected by 
their high equivalent width also comprised objects with 
very faint continuua, that would have fallen below the 
magnitude threshold of current Lyman break surveys, and 
also included two objects that could not be detected in 
Keck imaging of the optical continuum (1 a B = 27.8, 
V = 27.5, and I = 25.8; W\ > 400 A). 

In the present Letter we first present spectroscopic fol- 
lowups for the narrowband candidates of Paper I. In a 
small fraction of the cases, the spectrum shows both Lya 
and Civ A 1550 A, confirming the redshift identification 
but suggesting AGN-like properties. However, the major- 
ity of the spectra show only a single strong emission line. 
The absence of other detectable features, in combination 
with the high equivalent width of the selected candidates, 
identifes the single line as redshifted Lya emission, and 
argues that the equivalent width criterion (W\ 3> 100 
A) can, in fact, be used as a good diagnostic of high-z 
Lya-emitting galaxies. We then (§3) present results of 
a second deep narrowband search with a 78 A bandpass 
A ~ 6741 A filter (Lya at z ~ 4.54) and followup spec- 
troscopy, that confirmed two Lya-emitting galaxies at this 
higher redshift. In §4, we describe a very deep blank- 
field spectroscopic search (6 hr LRIS integration covering 
AA ~ 5000 - 10000 A) that yielded four emitters at red- 
shifts 3.04 - 5.64. Finally (§5), the data on emission- line 
objects from the imaging surveys in the two redshift inter- 
vals are combined with various spectroscopic surveys at 
lower redshift to show the evolution of the emission-line 
fluxes with redshift. We argue that the Lya-emitting ob- 
jects are significant contributors to the integrated star for- 
mation of the galaxy population throughout the z = 3 — 6 
redshift range, and that the integrated star formation 
rate of the Lya selected objects is flat, or possibly ris- 
ing through this redshift range, with a value greater than 
0.01 Mq Mpc" 3 for qp = 0.5 and H = 65 km s" 1 Mpc^ 1 . 
At the highest redshifts, most of the star formation may 
be occurring in objects of this class. 



2. SPECTROSCOPIC CONFIRMATION OF CANDIDATE 
LYMAN ALPHA EMITTERS 

Paper I provides a sample of 12 objects in the SSA22 and 
HDF fields, with W x > 100 A in the 5390/77 A filter. A 



further 7 objects were detected in a similar field surround- 
ing the z ~ 4.53 quasar BR0019-T522 (Storrie-Lombardi 
et al. 1996), giving a total sample of 19 objects over an 
area of 75Lj'. All the SSA22 and HDF candidate objects 
have now been spectroscopically observed using LRIS on 
the Keck II telescope, as have three of the objects in the 
BR0019-1522 field. Multi-slit masks with l'/4 slits were 
used, giving a resolution of 16 A (Cowie et al. 1996), and 
exposure times ranged from 1—3 hrs, with most objects 
having the longer exposure. Two objects in the HDF field 
have strong Civ emission, placing them at z = 3.40 but 
suggesting AGN activity (Fig. ||). These objects appear 
similar to the brighter objects in Pascarelle et al. (1996). 
However, the remaining objects have only one strong emis- 
sion line, and the continuua are too weak for absorption 
line identifications. 
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Figure 1. Spectra of strong Lya emitters in SSA22, LAI and LA3, (top 
two panels) compared with the spectrum of an [O II] emitter with ex- 
tremely high equivalent width (third panel) shifted so the [On] emission 
lies at the narrowband wavelength (5390 A). We also show for compar- 
ison one of the two Lya emitters (LAI) in the HDF field that shows 
strong Civ emission (bottom panel). The observed equivalent widths of 
the two SSA22 objects arc 325 A for the topmost object and > 845 A 
for the second. The observed equivalent width of the [O II] line in the 
comparison object is 220 A. The H/3+[Oni] complex is clearly evident in 
the case of the [O n] emitter, and these cases may be easily distinguished 
from genuine high-z Lya emission objects. 

A primary concern in the simple equivalent width selec- 
tion procedure is whether we are misinterpreting emission 
lines such as [O n] 3727 A and Ha 6563 A, that are common 
features in galaxies at moderate redshift, as higher redshift 
Lya emission. For objects where the continuum is strongly 
detected it is possible to use color information to diagnose 
the different lines (Paper I), but in cases in which only 
one emission line is seen, we can only make the distinc- 
tion based on the equivalent width strength of the line, 
assuming that an observed equivalent width W\ » 100 
A is generally Lya rather than unusually strong [O n] or 
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Ha. Spectroscopic followup is therefore essential to decide 
whether this procedure is reliable. 

It is easiest to argue that single-line objects are Lya 
emitters at lower redshifts, where the number of possible 
contaminating lines is smaller. An emission line at 5390 
A lies below rest-frame Ha, and only slightly above the 
[OlIl]+H/3 complex at 5000 A, so that the only probable 
sources of confusion are [On] 3727 emitters at z ~ 0.45. 
However, such objects would then have H/3 4861 and [O in] 
5007 lying at 7030 A and 7241 A, and we would ex- 
pect these emission lines to be strong, given the necessary 
strength of the [O n] 3727 line. In Figure |l| we compare the 
spectra of two strong emitters (W\ — 325 and > 845 A) 
found in the SSA22 field with one of the strongest equiv- 
alent width [Oil] emitters (W\ — 220 A) found in the 
Hawaii Deep Field Survey of Cowie et al. (1996), and 
with one of the objects with Lya and C iv emission de- 
tected in the HDF. The [O n] comparison spectrum has 
been wavelength shifted so that its [O n] line lies over the 
emission lines in the objects. It can be clearly seen that the 
absence of the H/3+ [O in] complex in the SSA22 spectra 
rules out the possibility that the objects are [O n] emitters, 
and leaves us only with the Lya interpretation. 
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Figure 4. The plots show Lya emission-line objects recovered from deep 
multi-slit spectroscopy of a field around the quasar BR1202-0725 using 
LRIS. The 4 detected objects have observed equivalent widths of > 600, 
> 600, 450, and 350 A. See Fig. 3 for the corresponding two-dimensional 
spectra. If the emission line were [O II], the H/3-t-[0 III] complex would lie 
in the 1600 - 1630 A region. 

At higher redshifts, the most problematic contaminants 
are likely to be high ionization extragalactic Hn regions, 
with strong lines of [O in] and Ha. Spectroscopic followups 
of candidate emitters from narrowband surveys where the 
filter bandpass lies in regions devoid of strong night sky 
lines, as is the case for the 5390 A and 6741 A filters, 
can rule out [Oin] by the doublet line structure. Ha 
from high ionization systems can have observed equiva- 
lent widths of several hundred angstroms (cf., Stockton & 
Ridgway 1998), and requires broad wavelength coverage of 



the [Oni] lines and/or the measurement of a continuum 
break across the emission line to confirm it as Lya. We 
do not find such contaminants evident among the narrow- 
band and spectroscopic samples discussed in the present 
paper. 



3. LYMAN ALPHA EMITTERS AT Z = 4.52 

Narrowband observations of the SSA22 field were ob- 
tained on the nights of UT 1997 September 21-22 with 
the LRIS camera on Keck II, through a narrowband in- 
terference filter centered at 6741 A, with a bandpass of 78 
A and transmission of 81%. The net exposure of 2.7 hrs 
was obtained as a sequence of 900 s exposures taken with 
10" offsets, with FWHMs of 0'.'85 (1.2 hrs) and 0"65 (1.5 
hrs). A corresponding 1.2 hr i?-band image was obtained 
as a series of 360 s exposures to define the local contin- 
uum. Three W\ > 100 A objects were detected, which are 
marked on the narrowband image in Figure El (Plate 1), 
and shown as 2-dimensional spectra in Figure [| (Plate 2). 

At this narrowband wavelength (6741 A) we may be 
contaminated by [On] 3727 at z = 0.8, Ha at z = 0.03, 
or [Oiii] or H/3 at z ~ 0.37. However, our spectroscopic 
followup again easily distinguishes such cases. The low- 
est flux object, with the weakest equivalent width con- 
straint (W\ = 120 A), is identified as an [On] emitter 
at z — 0.814, but the stronger object with substantial 
red wavelength coverage does not have detected [O in] and 
does not permit this interpretation. For the second strong 
object there is a strong break across the emission line, 
consistent with the Lya interpretation. Nor does either of 
the objects have emission lines near the primary line, that 
would permit the objects to be identified as [O in] or H/3. 
Once again there appears no choice but to interpret these 
lines as Lya. 
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Figure 5. Composite plot of Ha flux times rcdshift squared versus red- 
shift. At low rcdshift we show the /(Ha+[N II]) X z 2 fluxes (filled squares) 
from the z < 0.045 sample of Gallcgo ct al. (1995), converted to corre- 
sponding He* flux by multiplying by 0.75. The crosses show normalized 
Ha fluxes or converted, normalized [O II] fluxes from the Hawaii wide- field 
surveys (Songaila et al. 1994). The filled triangles show renormalizcd, 
converted [On] fluxes from a B — 25 spectroscopic sample (cf., Cowie 
ct al. 1997) with open squares showing renormalizcd /(Ha+[N n]) X z 2 or 
with circles showing renormalizcd /(Lya) at the highest redshifts if the 
spectrum docs not cover [O ll]. Finally, the remaining circles (at z ~ 3.43 
and z ~ 4.54 show /(Lya) X 2 2 /8.7 for the samples of objects detected in 
the two narrowband filters. The solid (qq — 0.5) and dashed (qo — 0.02) 
reference lines show the fluxes corresponding to star formation rates M 
= 10, 1, and 0.1 Mq yr" 1 for H Q = 65 km s -1 Mpc -1 and the M 
normalization discussed in the text. 
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The surface density of emitters from the narrowband 
search is sufficiently high (~ 15, 000/ 0° /unit z at z = 3.4, 
based on the 3 fields) that searches of existing very deep 
long-slit or multi-slit LRIS observations should yield such 
objects over the optical redshift range (z ~ 2 — 7). In- 
deed, Franx et al. 1997 identified a second z ~ 4.92 object 
by inferring that the single strong emission line seen in a 
serendipitously placed slit during spectrosopic exposures 
of their color-selected, gravitationally lensed galaxy, was 
also due to Lya at their target object's redshift. In our 
own archival LRIS data the best suited exposure is that 
around the quasar BR1202-0725, where a 6 hr integra- 
tion in a single mask configuration was obtained under 
photometric conditions and with reasonable seeing (UT 
31 March - 2 April 1997; 4 hrs obtained with 0'.'8 - Of! 85 
FWHM, and the remaining 2 hrs with - 0"95 FWHM). 
The data were taken to study color-selected and emission- 
line-selected objects in the quasar field (Hu, McMahon, & 
Egami 1996; Giallongo et al. 1998) and the mask com- 
prised only twelve slits, with the blank sky area cover- 
ing roughly 200 arcsec 2 . A search of this spectral im- 
age yielded four extremely high equivalent width emitters 
ranging in redshift from z — 3.04 — 5.64. These spectra 
are shown in Fig. |^ and in Plate 2. The two highest red- 
shift objects can be seen in the red continuum, and have a 
strong break across the emission line, consistent with the 
expected Lya forest decrement. 

The surface density of emitters in the spectroscopic field 
is roughly 60,000 deg~ 2 /unit z in the accessible redshift 
range, consistent with the narrowband searches when ac- 
count is made for the slightly fainter flux limit in the spec- 
trocopic data (1 er <~ 10~ 18 ergs cm -2 s _1 ) and the current 
uncertainty in the numbers. Ultimately, the census of the 
populations is best done with the narrowband data with its 
uniform selection procedures and well-defined flux limits, 
but the present data suggest that we will see substantial 
populations of these objects to very high redshift. 

5. DISCUSSION 

Since resonant scattering enhances the effects of extinc- 
tion, it is harder to convert the Lya emission into a massive 
star formation rate than it is for line luminosity diagnostics 
such as Ha and [On] 3727. For the present calculation, 
we assume that extinction may be neglected in comput- 
ing the required massive star formation rates, which then 
constitute a minimum estimate. However, upward correc- 
tions to this value are unlikely to be larger than a factor 
of two, since the observed rest-frame equivalent widths lie 
in the 100 — 200 A range — close to the maximum val- 
ues that are obtainable from ionization by a massive star 
population (Chariot & Fall 1993). Then, assuming case 
B recombination, we have L(Lya) = 8.7L(Ha) (Brock- 
lehurst 1971), which using Kennicutt's (1983) translation 
of M from Ha luminosity, gives M — (L(Lya)/10 42 ergs 
s _1 ) Mq yr" 1 . In order to cross-calibrate to [On] fluxes 
at lower redshift we assume /(Ha+[N n]) = 1.25 /([On]) 



based on the mean values of the ratio in both the Gal- 
lego et al. (1995) and Hawaii Deep Survey (Cowie et al. 
1996) samples. We also assume /(Ha+[N nj) = 1.33 /(Ha) 
(Kennicutt 1983). 

In Figure ^ we compare the range of Lya luminosities to 
the range of line luminosities in lower redshift objects. The 
plot shows the quantity z 2 f vs redshift, where Ha+[Nn], 
[O n] , and Lya fluxes have been converted to Ha fluxes 
using the relationships above, and we have restricted our- 
selves to the imaging data in which the fluxes of the Lya 
are well determined. The comparison objects are drawn 
from the surveys of Gallego et al. 1995 (filled boxes), 
Songaila et al. 1994 (pluses), and the Hawaii B — 25 sam- 
ple [open boxes (Ha+[N n]), triangles ([On]), and circles 
(Lya)]. The solid (q = 0.5) and dashed (qo — 0.02) lines 
on Fig. H show the fluxes corresponding to stellar mass 
production rates of 10 Mq yr^ 1 , 1 Mq yr -1 , and 0.1 M Q 
yr^ 1 for Hq = 65 km s _1 Mpc" 1 . Maximum formation 
rates locally are around a few Mq yr -1 , rising to values 
of just over 10 M© yr -1 above z = 0.6. The Lya fluxes 
at the higher redshifts are then consistent with this value 
or just slightly smaller depending on the extinction cor- 
rection. The CADIS results (Thommes et al. 1998) would 
lie a factor of several times higher in flux than the two 
Hawaii filter samples, but both Keck spectroscopy and re- 
peat Fabry-Perot observations have disproved the original 
z ~ 5.7 candidate selection (Meisenheimer 1998). 

The minimum integrated star formation rates at z = 3.4 
and z = 4.5 are 0.006 Mq Mpc -3 yr" 1 and 0.01 M Q 
Mpc~ 3 yr~ x respectively for go = 0.5 where the first value 
is slightly smaller than that given in Paper I since AGN- 
like objects are excluded. Both values are lower limits 
calculated in the absence of extinction and the z — 4.5 
value is based on a single field. Within the substantial un- 
certainties of the as yet small samples, the results suggest 
that the star formation rates in the strong emission line 
population are constant or may possibly be increasing with 
redshift from z = 3 — 6, in contrast to color-based samples 
where the rate is declining at higher redshifts (Madau et 
al. 1996, 1998). This is consistent with the broad ex- 
pectation that as we move to higher redshifts and earlier 
stages of galaxy formation, an increasingly larger fraction 
of the star formation should be in strong Lya emitters that 
correspond to the youngest galaxies. 
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Fig. 2. — (Plate 1) A 2.7 hr Keck narrowband image of the 6' x 4' region around the SSA22 field through a 78 A filter centered at 6741 A. Objects with observed 
equivalent widths W\ > 100 A are circled. The object with the lowest observed equivalent width (W\ = 120 A) is spectroscopically identified as a z = 0.813 
[O n] emitter. However, the two higher equivalent width objects appear to be Lya emitters (W x = 200 A, z = 4.521 and W\ = 1000 A, z = 4.528). 



Fig. 3. — (Plate 2) Gray scale images of emission-line objects selected in the spectroscopic observation of the BR1202-0725 field and of the two emitters in the 
SSA22 field (Plate 1) selected in the 6741/78 A filter (top two panels). The images have been registered in wavelength with a range from 4400 A to 9700 A. The 
lower seven panels show emission-line objects at z < 1.20 in the spectral image. For each object, the [On] A 3727 line is indicated. The two middle panels show 
the z = 4.69 quasar, BR1202-0725, observed at two different orientations. Narrow, spatially separated Lya emission without strong continuum is seen at both 
angles, corresponding to the two companions to the quasar (Hu, McMahon, & Egami 1996; Fontana et al. 1996; Petitjean et al. 1996; Hu, McMahon, & Egami 
1997). The next four panels show the single emission-line objects (thick tick marks) in the BR1202-0725 field. The two longer wavelength objects at z = 4.19 
and z = 5.64 have continuum breaks across the emission line of 2.3{| q and 2.6{f| respectively, suggesting that this line is Lya, while the two lower wavelength 
lines, where no continuum is visible, have none of the additional emission lines which would be seen if they were at lower redshift. (The longest wavelength 
emission line object lies close to a bright foreground galaxy at z = 0.236, whose spectrum can be seen at the bottom edge of the image. However, the required 
velocity separation of —1500 km s _1 would seem large for the line to be interpreted as Ha from an extragalactic Hu region associated with this object.) For 
the SSA22 (z ~ 4.52) objects, the spectrum in the upper panel shows a continuum break, while the lower panel's spectrum is only an isolated emission line. 



